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Abstract. The velocity of a passive particle in a one-dimensional wave 
field is shown to converge in law to a Wiener process, in the limit of 
a dense wave spectrum with independent complex amplitudes, where 
the random phases distribution is invariant modulo %j2 and the power 
. spectrum expectation is uniform. The proof provides a full probabilis- 

tic foundation to the quasilinear approximation in this limit. The result 
extends to an arbitrary number of particles, founding the use of the en- 
semble picture for their behaviour in a single realization of the stochastic 
wave field. 
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1. Introduction 

We first recall the physical setting in sec. 11.11 Indeed the quasilinear ap- 
proximation is an ubiquitous scheme for deriving irreversible, diffusion-like 
equations from many-body dynamics, involving a "propagation of chaos" 
kind of argument in a system with mean-field behaviour. Original deriva- 
tions of this approximation are sketched in sec. 11.21 A different, more re- 
cent analysis in the framework of hamiltonian chaotic dynamics is recalled 
in sec. 11.31 As these arguments are well detailed in the literature, we do not 
reproduce the calculations and proofs but merely highlight their key points. 
Our main result and particulars of the present work are stressed in sec. 11.41 

1.1. Physical setting. The motion of a particle in the field of many waves 
is a fundamental process in collisionless plasma physics. Even if the particle 
motion does not feed back on the wave parameters, viz. for a test particle, 
undergoing passive transport, this problem still presents open issues. Its 
most elementary instance, in one space dimension, is also a benchmark for 
approximation techniques. 

This one-dimensional problem describes the motion of a particle in a lon- 
gitudinal, electrostatic, time-dependent potential. Electrostatic modes occur 
in various contexts |GoRu95, DMA05J, including (i) the non-relativistic 
regime of Coulomb plasmas, where magnetic fields are negligible ; (ii) par- 
ticle motion parallel to the applied magnetic field in strongly magnetized 
plasmas ; (iii) particle motion along the axis of a waveguide, such as travel- 
ing wave tubes used as amplifiers in telecommunications. The time depen- 
dence of the field leads to the propagation of waves, which are longitudi- 
nal : Langmuir waves are the simplest collective modes in plasmas. When 
it applies (in particular for hot plasmas), the neglect of collisions in particle 
dynamics within plasmas rests on the long-range nature of Coulomb inter- 
action leading to a mean-field picture in both the Vlasov kinetic equation 
and the Euler fluid models. 

In many situations, the wave field evolution involves a response to the 
particle motion. However, in some instances the particle feedback on the 
electrostatic field is negligible, and one may take the field as given. The 
equations of motion for a particle with charge c and mass m then read 

§ - m 

(2) f = 1 *(*(.),,) 

where the electric field E is a prescribed process. It is convenient to rep- 
resent E as a Fourier series, E{x,t) = Y, m E m e l ( kmX ~ a>mt+(p ' n \ where ampli- 
tudes E m and phases (p m may be tunable, while (k m , G) m ) are given by the 
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waves dispersion relation^ The key effect of a single wave with phase ve- 
locity V(p : , n = COm/km on a particle with velocity v is a tendency [DEM05J 
to reduce the relative velocity |v — v<p. m |, and this effect works best when 
ttl(v — Vm , m ) 2 ~ \zE m /k m \. The competition between two waves m,m' in 
attempting this synchronization is measured by the wave overlap parameter 

e I 1/2 
m 

When this parameter is small, a particle cannot interact strongly with both 
waves simultaneously, and the dynamics can be analyzed perturbatively. 
Actually, the dynamics ®-0 is we U known to be nonintegrable as soon as 
there is more than a single wave phase velocity ; the two-wave model is a 
paradigm of hamiltonian chaos with 1.5 degrees of freedom, with a KAM 
limit as s — > 0, and transition to "large scale chaos" as s > 1 . r Es 8 51 IeTEsQ3H 
Denoting by Av<p the typical relative velocity of a wave with respect to 
its nearest neighbours, by k typ a typical wavevector, and by E typ a typical 
amplitude, the regime of interest for this paper is the dense spectrum limit, 
where a particle is typically influenced significantly by many waves ; in this 



regime the typical resonance overlap parameter s = 4y | (cE tyv ) / (mk tyv )\/ Avq, 
is large. It is then tempting to consider the acceleration in the right hand 
side of © as an approximate white noise, and the particle velocity V as a 
kind of diffusion process : this is the core of the quasilinear approximation 
HRoFi61l IVVS611 IVVS621 IDrPi62L The latter is widely used, in diverse 
physical contexts, as it is easily implemented and relies on simple ideas, 
which we comment in the following. 

1.2. Original derivations of the quasilinear approximation. Classical 
derivations of the quasilinear approximation in plasma physics textbooks, 
e.g. HKad65[|GoRu95 , HaWa04], start from viewing the motion of the test 
particle as the transport of a measure d/i = fdxdv on (x,v) space (with 
/(.,.,?) possibly a distribution), 

(4) d t f + vd x f=--E(x,t)d v f, 

m 

and begin an iterative solution with respect to E, 

e f f 

f(x,v,t) = f(x-vt,v,0) / E(x-v(t-t 2 ),t2)d v f((x-v(t-t2),v,t2)dt2 

t t m 

+(-) 2 / I 2 E(x-v(t-t 2 ),t 2 )E(x-v(t2-ti),t l ) 
m Jo Jo 

(5) d v f((x-v(t-t l ),v,ti)dtidt 2 . 



^Note that this model differs from stochastic acceleration problems in a random poten- 
tial, for which the field E(x,t) reduces to a static random E (x) . 
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On performing an x average which highlights the correlation function of 
the electric field E, one then relies on independence of E from the (slaved, 
passive, tracer) particle distribution d/i to eliminate the first order term, 
and obtains an integro-differential evolution equation for the jc-averaged /. 
Then, on considering that the velocity process V must be Markov on time 
scales longer than the correlation time of E, the equation for /(v, t) reduces 
to 



£>( v ) = -!_ / {E{x,t)E{x-vx,t-x))dx= 7 ^ T [ 8((D-kv)(\E k \ 2 )dk 



is determined by the wave field lagrangian autocorrelation, with appropriate 
averaging (•) and assuming that phases <p m are independent and uniformly 
distributed. The Fourier form in ©, with a Dirac distribution, obtains in 
the continuous spectrum limit. 

The "appropriate averaging" (•) may imply (see e.g. sec. 9.4 in [HaWa04J) 
that one no longer considers the evolution of test particles velocity distribu- 
tion / f(x, v, t) dx in a single realization of the field E but rather the expecta- 
tion of this J f(x, v, t) dx with respect to the ensemble of wave fields. Such a 
view pertains to the statistics of particle velocities collected from repeated 
experiments, but it does not apply a priori to the description of transport in a 
single realization, as stressed in more general terms e.g. p. 45 in [HaWa04J. 

This derivation may be criticized (within its own viewpoint) on the ground 
that, however small the coefficient E may be, the differential operator d v 
is unbounded for many function spaces. Formal, diagrammatic HBo62al 
IBo62b[ [Bo65 , ThBe73 , BrFr74J expansions in E d v are therefore less straight- 
forward than they may seem. 

An alternative derivation, based on particle motion and E -power expan- 
sion, also leads to the diffusion equation © via its Langevin counterpart, 
assuming that the particle velocity is a Markov process and computing the 
first two moments of its increments [St66J. In this context, the "random 
phase approximation" is actually invoked so that, to practical ends, "la- 
grangian" (as seen by a test particle) phases k m (x — v(t — tj)) + (p m — G) m tj, 
can be considered as independent (uniformly distributed modulo lit) ran- 
dom variables for any relevant sequence of times tj and wave indices m, 
viz. not only at a single time (this is imposed by the very distribution of 
parameters <p m ) but as if their values were "refreshed" repeatedly. The bold- 
ness of such an assumption, akin to the propagation of molecular chaos in 
gas theory HKac56[ |Kac591 , fueled the debate on the validity of the quasi- 
linear approximation (as a preamble to the further debate focusing on the 
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self-consistent problem, where wave amplitudes and phases evolve under 
particle feedback) HCEV901lLXW93llLaPe99ll . 

Mathematically, the "repeated random phase approximation" is valid, un- 
der a few more technical conditions [Pa Ko74H . in the limit e — \ after a 
time rescaling, x = £ 2 t, when the field E is mixing, in the sense that the 
process E(.,t) is adapted to a family of a-algebras < s < t < oo, with 
&s\ ^= f° r ^^S2<si<ti<t2<°°, with a probability measure P such 
that the rate function 

(8) p(?):=sup sup |P(A|fl)-P(A)| 

s>QAe<0?£_ t ,Be3?% 



satisfies the condition JJ° Wp(t)dt < <*>. Typical examples of such mixing 
processes £ are ergodic Markov processes on a compact space [PaKo74|, 
but time-periodic fields as discussed e.g. in Refs HCEV90i rBeEs98a, ElEs03 , 
ElPalO] fail to meet the mixing condition. 

The use of an adjoint formulation instead of trajectories is generally mo- 
tivated by the traditional viewpoint of kinetic theory, interested in following 
many particles (in which case, including the self-consistent dynamics where 
the evolution of E depends on /, measures provide a natural description, see 



e.g. ch. 1.5 in [Sp91 1), by the fact that the Vlasov and diffusion equations 
are linear for /, and by the familiar description of Markov processes in 
terms of their generator. Yet a single physical realization of the wave field 
E acts on a particle distribution quite differently from the way an ensem- 
ble of independent realizations would act on a single particle [BeEs98b |. 
The decorrelation assumption is crucial in claiming that the ensemble may 
describe a single experiment with many particles. Besides, if the Markov 
assumption fails, the single-time distribution function f(x, v, t |jco, vo, to) may 
fail to describe properly the joint n-time distribution F(xi,vi,t\ . . .x n ,v n ,t n ). 
Therefore we revisit the derivation of quasilinear equations from a particle 
viewpoint, and possibly reach a Markov description in an appropriate limit. 

1.3. Hamiltonian dynamics approach. This dynamics-based program was 
significantly advanced by Benisti and Escande [BeEs97, BeEs98a|, who 
proved the validity of the velocity diffusion picture for the dynamics de- 
fined by hamiltonian 



(9) H=^—+srf cos(q - mt + (p m ) 



M 



2m 



-M 



in the limit M 3 / 2 3> /m — > °°, when phases (p m are independent and uni- 
formly distributed in [0, lit] . Their derivation relies on the strong chaos (as 
s — > oo) in particle dynamics associated with the limit, and on the fact that, at 
a time?, only waves with a phase velocity such that |v<p — p(t)/m\ < Avb act 
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strongly (nonperturbatively) on the particle. Waves beyond the "resonance 
box half-width" Avb ~ 5(«sz/ /m) 2//3 can be eliminated from the dynamics 
(their overall statistical effect is exponentially small in |v<p — p/m|/Avb) by 
a canonical transformation, so that the velocity process is Markov on scales 
wider than the resonance box. On the other hand, for shorter time scales, 
the particle velocity needs a time of the order of unity to sample correla- 
tions associated with the discreteness of the frequency spectrum, so that it 
is chaotic and wanders so much that it eventually moves to another reso- 
nance box. Moreover, for short time scales, they show how to relax the 
assumption that all phases are independent to the requirement that any two 
phases influence negligibly the particle motion [BeEs97 , BeEs98a, ElEs03 1 . 

This argument was complemented by the observation that the short-time 
quasilinear approximation holds for times < t < D~ 1//3 lns, and that the 
Markov approximation holds for times t > D l / 3 HEsE1021 lEffisOl (IT 1 / 3 
is also related to the Lyapunov time scale for the divergence of microscopic 
trajectories in a typical wave field [ElEs03J), so that the quasilinear approx- 
imation holds for all times in the dense spectrum limit s — > °°. 

For technical simplicity, the hamiltonian model © involves three re- 
strictions with respect to the original dynamics ©-© : all amplitudes are 
equal, all wavevectors are equal, and all phase velocities are equally spaced. 
Benisti and Escande [BeEs98a| sketch how their arguments can be extended 
to lift these restrictions. The hamiltonian © also stresses the spectrum dis- 
creteness time scale, as Av<p = 1, which can generate correlations over long 
times HBeEs971lBeEi98bl . 

1 .4. Position of this work. In the present work we extend the approach 
initiated in [E107J, lElPalOM and revisit the Benisti-Escande result with the 
language of probability theory. We express the wave field as a sum of /V — > 
oo independent components per unit frequency interval, so that the overlap 
parameter s diverges in the limit /V — > °°. We first prove in Theorem l3.4l that. 
in the resulting dense wave spectrum limit s — > °°, the wave field acting on 
a particle for < t < lit converges in law to the field associated with a 
"white noise". This enables us to derive Proposition 14.21 and Theorem 14.31 
on particle motion. 

Proposition 14.21 shows that, for M — » °°, for fixed wave power spectral 
density with /V — > °° so that s — >■ oo, the velocity of a single particle in 
the wave field converges in law to a Wiener process over the time inter- 
val [0, 111] . While Benisti and Escande emphasize a hamiltonian dynamical 
system approach, we focus on the velocity and express our limit theorem as 
a convergence in distribution result, following essentially from central limit 
averages on the wave field. The convergence in distribution was clearly 
understood in [BeEs97, BeEs98a|, in particular through the statement that 
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the influence of waves outside a resonance box is only perturbative on the 
statistical properties of the dynamics (p. 914 in [BeE s98alO . The focus on v 
is also central to the arguments in [ElEs03J which involve the characteristic 
function $(7,?) := E exp(iy(v(f) - v(0))). 

We also pay attention to the behaviour of an arbitrary number jV of par- 
ticles moving in the same wave field. Their evolutions are not independent 
processes for finite srf , so that the diffusion equation © does not describe 
the evolution of the empirical distribution ^^^Y^iLi <5(v — pi/vcv) ; this was 
stressed in [BeEs98bJ. However, our main result, Theorem l4.3l proves that, 
in the limit s — )■ °°, particles do diffuse independently, even in the same 
wave fielcfl - in agreement with the view that they generally are in different 
resonance boxes. Thereby we extend to a broad class of wave fields the 
conclusion of [ElPalOJ, which assumed a wave field generated by Wiener 
processes (viz. the fields obtained in the dense spectrum limit). This result 
provides some support to the traditional view that a single realization may, 
in some cases, be approximated by an ensemble. 

This paper is organized as follows. We state our assumptions on the 
wave field in section [2l These enable a fast proof of the convergence of 
elementary processes associated with the wave field in section [3] Thanks 
to Ref. [ElPalOJ and the continuous mapping theorem [KalOl], the con- 
vergence of the ,/K-particle velocity process to the diffusion limit follows 
immediately in section |4l We stress the interpretation of our techniques and 
results in section[5l We conclude with a discussion of open issues. 

2. Wave field assumptions 

A random variable (r.v.) a is symmetric [Kah85, KalOlJ if a and —a 
have identical distributions. Then Ea* = for odd k if the expectation 
exists. 

ASSUMPTIONS 2.1 (S2, S4). Given M G N = {0, 1,2,...} and N G N = 
{1,2, . . .}, consider (2M+ \)N complex random variables a m ,n =A mjn e i<Pm '". 
We say that the (X m , n 's meet assumptions (S2) if 

(1) they are independent and symmetric, 

(2) EAl n = I, 

(3) sup OTn EA^ n < C4 for some C4 > 1. 

We say that they meet assumptions (S4) if, moreover, the r.v. a] nn is also 
symmetric. 



The contrast between this conclusion and Benisti-Escande's |BeEs98b| might be at- 
tributed to the asymptotic nature of our result, as s — > 00, We do not provide estimates for 
the "convergence rate" of the empirical distribution to its Fokker-Planck limit. 
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The additional condition for (S4) may be called "four- symmetry" for the 
r.v. a. Examples are (i) the r.v. e 1 (c+ Kn / 2 ) t with fixed c and F(K = k) = l/4 
for k E {1, 2,3,4}, (ii) an isotropic complex r.v., viz. a =Ae i<p such that <p is 
uniform on [0, 2it\ (which corresponds to a Steinhaus sequence (p m . n / (2n) 
[Kah85]) and independent from A, and in particular (iii) a standard normal 
complex r.v. (isotropic, with exponentially distributed A 2 ). 

REMARK 2.2 . The third condition in (S2) is unnecessarily stringent ( though 
being met for many physical cases), and could be relaxed to a Lindeberg- 
type condition. 

Occasionally we identify M 2 with C to minimize the amount of notations. 
We denote by B the standard brownian motion in C(M + ,IR) and by W the 
standard brownian motion in C(M + , C), so that B, y/l^iW and \/2ZW are 
independent and identically distributed (i.i.d.). 

3. Convergence of the controlling wave processes 

Given N real parameters a n E [0, 1] (1 < n < N), we first introduce the N 
complex- valued processes, for 1 < n < N, 

1 r* M 

(10) u M {t)= £ a in . n e-« m+ ^ s ds 

V2n Jo m f- M 

for t E R. By construction, u^ is analytic for any finite n,M, and e lfV duff / dt 
is a family (1 < n < N) of independent 27T -periodic complex processes. In 
the limit M — > °°, the processes u^ lose their smoothness (as, typically, their 
Fourier coefficients decay slowly), but Proposition 13.21 shows that they al- 
most surely (a.s.) admit a Holder continuous limit u n . 

Specifically, we characterize the smoothness of a function y E C(R, C) 
by its modulus of continuity [Kah85, KalOTTl. 

(11) tty:]0,oo[-»[0,oo] CDy(h)= sup \y(t') -y(t)\ . 

\t-t'\<h 

Our first objective is a gaussian convergence theorem, in the limit Af — > °°, 
for the complex- valued process = N~ l l 2 Y<n u n- Let 

02) ^(o=4=ix« 

for ?Gl. Note that if the a„'s do not vanish and the a„ v , are i.i.d., processes 
u n are not i.i.d., but they remain independent with closely related moments. 
Forge C^QO^j^C), let 

f7.lt 

(13) 8m,n = {^r 112 g^C-'^+^'dt. 

Jo 
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We also introduce the N complex- valued processes, for 1 < n < N, 

M 



1 f' M 

(14) yM(t) = — £ a nhn z-» ns ds 

Jin Jo 



m=—M 

M 



for t e M. In case the a m ^ n are i.i.d., the processes y„ are i.i.d. for given M. 

PROPOSITION 3.1. Let (768. Ify e C(M,C) has a modulus of continuity 
COy, and u E C(R, C) is defined by u{t) = Jq e~ ias dy(s) for (el, then its 
modulus of continuity satisfies CO u (h) < (1 + \o\h)(Oy{h) for h > 0. 

Proof First note that, for any t,t' e M, 

e im \u{t')-u(t)) = J e- ia ^dy(s) 

= y( t ')-y(t) + -l)dy(s) 

= y(t') -y(t) + [(e- 1 ^') - l)(y(,) -y(f))' 
(15) + ^ -y(0)iae- lcJ ( iW ')d^ 



j t 



by integration by parts. The middle term in the right hand side of (1151) 
vanishes, and we estimate the sum using triangle inequality for t < t', 
(16) 

\u(t')-u{t)\ = \^\u{t')-u{t))\ < W) -y{t)\+f t \y{s) -y{t)\\o\ds 

from which the claim follows by definition of the moduli of continuity. □ 
For < /3 < 1 and p e N, we denote by C P ^(R, C) the class of continu- 
ous complex-valued functions of a real variable, with p continuous deriva- 
tives, such that their p-th derivative is Holder continuous with exponent /3. 

PROPOSITION 3.2. Let u„ and yf? be defined by $dU$ and (O under as- 
sumptions (S2). For any < j8 < 1/2, and for any n, the sequences y% and 
uff converge a.s. in C ^ (R, C) as M — > °°. 

Proof The y statement results immediately from Theorem 3, Sec. 7.4 in 

■ . r, o/'+ii — ? ? r -i-i 2^ +1 — 1/2 

HKah85L as we compute the sums s ■ = )'~ ,.- m Z EA^„« [m J 2/— 1/2 ? 
2-7+i/2 f us i n g the fact that EA^ n = 1. 

The u statement then follows from Proposition l3.1[ □ 

PROPOSITION 3.3. Let u„ and y^ be defined by rfTDl) and (O under as- 
sumptions (S2). For any g £ C ([0,2^],M), consider the complex random 
variables {g,u%) := fg* g(t)dfl(t). 
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(i)ForanyM, E(g,u* f )=0, E(g,u^) 2 = l^ = _ M g 2 m ^ n Ea 2 ^ n andE\(g,u^)\ 2 ■ 

Moreover, sup WjM E | (g, f < (2 + Q)||g||J. 
(//J Assume further that the a nijn 's are four-symmetric. Then as M —¥ oo, 
f/ie complex r.v.'s (g,u%) converge a.s. to a r.v. (g,u n ) such that E(g,u n ) = 
0, E( gl u n ) 2 = 0, E\( gl u n )\ 2 = fo n g 2 (t)dt, and sup ;7 E\(g,u n )\* < (2 + 

Q)IHll. 

Proof Calculations are straightforward as the given test function g is con- 
tinuous and [0, 2n] is compact : 



1 f 271 M , \ 

(17)E(g, M f) = — / g(t) £ Ea ra ,„e- 1 ( ra+(7 »"dr = 0, 

V 27t JO m =_ M 
1 /-2s r 2K M 

ns^f = / 8(t)g(s) £ Ea^ e-^)(^)d/d5 

2tt Jo Jo m f- M 



M 



(18) = £ E<„£„ 

m=— M 

r-2;r /-27T M 



- 1 /-2S /-27T m 

E|(g,«f)| 2 = X-/ / g(f)*W £ E^e-N^ 
2tt./o Jo 

M 

(19) = £ EA^^I 



Given that EA 2 n n = 1, the latter expression yields^by Parseval's identity 
(20) E|(g, M f)| 2 < £ |g m , n i 2 =/ s 2 (f)df=|Hl2 

jtj — — oo *■* 



with equality in the limit M — > <*>. Finally, 

= ^ £ a mi,nMm2,n&m2,n®m4,n8^,n8m 2 ,ng m3 jig m4j n 

2 



IS I 2 1 




\&mi,n\ 


8ni2,n | 



(21) < (C 4 + 2)||g|g ; 



As pointed out by a referee, this argument reduces to Bessel's inequality, when one 
views uff as a sum of 2M + 1 basis functions, in the Hilbert space (whose elements are 
stochastic processes u) with scalar product (u,v) = E J n u*(s)v(s)ds. Our assumptions 
on the r.v.'s a mj „ ensure orthonormality of our basis. 
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where the first equality follows from the definition of Fourier coefficients 
g m>n , the second equality from the known first two moments of a, and the 
final inequality from the bound C4 on EA 4 . □ 
Now we can prove our main claim, 

Theorem 3.4. Under assumption (S4), the process Uff defined by 4721) 
converges in distribution to the brownian motion in C( [0, 2n] , C) as N — > °° 
and M — )■ 00. 

Proof First, consider the process Un = lim^^ooL^f in C 0, ^(R,C) for any 
< j8 < 1/2. The convergence is a.s. since Uff is a finite linear combination 
of the processes u%. Given any g 6 C ([0,27r],R), we show below that the 
r.v. Zm '■= (g, Un) converges in distribution to a normal r.v. Z = X + iY with 
EZ = 0, EX 2 = EY 2 = ±\\g\\l, E(XY) = 0. AsC 1 is dense in L 2 , the same 
holds true for g e L 2 ( [0, lit] , R) , which will imply that the limit (g, U) is the 
Wiener integral [ Nu06ll . 

The first two moments of Z follow easily from the fact that the r.v.'s 
Cn = ^n + iln '■= (g,Un) are independent. Proposition 13 .3 1 states that E^ n — 
Er] n = 0, E (2^7]„) = 3 E (g, u n ) 2 = and E - T]„ 2 ) = 91 E (g, u n ) 2 = 0. 
Besides, E(<^ 2 + 7]„ 2 ) = E|(^, a „)| 2 = ||^|| 2 . 

The fourth moment condition implies the Lindeberg condition on the se- 
quence C, n (alternatively, one may adapt the standard proof of the central 
limit theorem via the characteristic function), so that N~ l l 2 Y^ = \ Cn con- 
verges in distribution to a normal complex random variable, by the gaussian 
convergence theorem (e.g. Theorem 5.12 in [KalCEQ). D 

Remark 3.5. Our statement holds for arbitrary choice of coefficients o n , 
essentially thanks to the fact that, for any o, functions (2^:) _1//2 e 1( - ra+c7 ) f 
form an orthonormal basis o/L 2 ([0,2^:],C). In the special case where the 
0C mjn are already normal, each u n is already a Wiener process. 

Remark 3.6. In the case where all o n = 0, the processes u n (t) — j^u n {2n) 
define In-periodic functions in C '^ (R, C) ; their restrictions to [0, 2%\ con- 
verge to the brownian bridge (see UKalOlll . ch. 13) and the N~ l l 2 Y*n a m,n 
converge to i.i.d. normal r.v.'s. 

Benisti and Escande [Ben95, BeEs97, BeEs98a| consider the case where 
CL m .n is uniformly distributed on the unit circle, and work with N = 1. Our 
statements do not formally apply to such a case. But they let their wave 
amplitude g/ — > °°, so that s — > °° and Avb — > 00 . To keep finite velocity and 
amplitude scales, we reformulate their case by relabeling the waves with 
integer- valued index m' = mN + n, letting o n = n/N, and rescaling time as 
t' = t/N, so that m't' = (m + o n )t. To address finite t' scales (of interest to 
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them), we need to extend our previous statements to arbitrarily large time ; 
the following statement is a first step in this direction. 

Proposition 3.7. Under assumption (S4), assume further that o n = n/N. 
Then the process defined by 4721) converges in distribution to the brow- 
nian motion in C(R, C) as N — >■ °° and M — > °°. 

Proof It suffices to prove convergence over an arbitrarily long time interval 
[0, 2nK], with K e No- To extend the previous theorem to K > 1, consider 
first a subsequence N = N'K with N ; — > °°. Then let s = Ks' and decompose 
n = n' + kN' with 1 < ri <N' and < k < K- 1. Note that 

(22) £7^(0 = r ! 4"± f^^e-^' 1 '^ 

n =\ v 2tt Jo m 

where the latter expression is equivalent to a process y/KU$ K (t /K), up to 
the K - 1 terms for which MK <mK + k< MK + K - 1. These if - 1 
terms do not spoil the limit as their contribution vanishes a.s. for N' — > °°, 
while by Theorem 13.41 the process y/KUj^, K (t / K) converges in distribution 
to y/KW(t/K), which is distributed as W(t). 

Now, if TV = N'K+k with 1 < k < K, the difference U$-{\ +k/N)-^ 2 U^, K 
converges a.s. to zero as — > °°, while 1^^/^00(1 +k/N)~ 1 / 2 = 1, so that 
the sequence Ujf converges like the subsequence U^, K . □ 



Remark 3.8. In rough paths terms (see e.g. [|Lej09[, sec. 8.4, and [FrVi06J 



for definitions and notations), Theorem \3.4\ corresponds to the natural ex- 
tension or lift Un (with Ujy = Un) converging in distribution to the geomet- 
ric enhanced brownian motion in C°'^([0,27r], G 2 (C))for 1/3 < /3 < 1/2. 



4. Particle motion 

We now turn to the solution of differential equations with control Ujf, 
viz. to the motion of particles in the wave field associated with the u„ 's. 
Since the latter functions are C , integration against them must be inter- 
preted so that the limit differentials d9t£7, d3£7 are the Stratonovich ones 
I WoZa65l IDo77[ ISu78l . This is satisfactory for the physicists applying e.g. 
diffusion models, but our result goes further : this formulation opens the 
way to analysing almost every single realization of the underlying noise, 
which need not be gaussian. 
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Specifically, the motion of a particle in the prescribed field of electrostatic 
waves is described by the system 

(24) < = -p^dt, 

m 

N M 

(25) dp% = AT 1/2 £ £ A m>n sin(^ (r) - (m + a n )t + (p m . n ) dt , 

n=\ m=—M 

(26) = sin(^(0)d9t<(0+cos(^(0)d3<(0, 

(27) 4(0) = g , P%(0)=Po, 

where srf is an overall amplitude scale (incorporating e) for the waves and 
m is the particle mass. We rescaled the particle momentum p by this overall 
amplitude to construct an appropriate limit below. 

Remark 4.1. In the special case where o n — for all n, the N -averaging 
generates gaussian coefficients for the Fourier wave components for each 
m. In the case where o n = n/N, the wave field has period 2nN, but its 
sampling over the shorter interval [0, 2k] prevents the observer in the limit 
N — > oofrom distinguishing it from an actual white noise. 

The previous section implies that, in the limit M — v °°,N — > °°, the equa- 
tions of motion may be interpreted as 

si 

(28) dQ = —Pdt, 

m 

(29) dP = sin(Q(f))od9l£/(0+cos(Q(0)od3£/(0, 

(30) Q(0) = q , P(0)=po, 

where od denotes the Stratonovich differential [Wo Za65l . 

Proposition 4.2. In the limit M -> o°,N ->■ °° the process (q^,p^), de- 
fined by rf2^-(|231)-ri2"7l) under assumption (S4) with (qo,Po) £ K 2 converges 
in law to (Q,P), where 

(31) GO) = q + — (p t+f t B(s)ds), 

m Jo 

(32) P(t) = p + B(t), 

with B the standard one -dimensional Wiener process in C( [0, 2k] , M). 

Proof Theorem 13 .41 ensures that the limit U is a standard complex Wiener 
process. Then, for the system d28T>- (|29l -(l30l) the mapping C([0,2rc],C) -> 
C([0, 27r],IR 2 ) : t/ i-> (P, 2) is continuous for the topology of uniform con- 
vergence HDo77llSu78L and the continuous mapping theorem [KalOlJ trans- 
fers the convergence in distribution from the control U to the particle evo- 
lution (P,Q). 
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The proof then follows Ref. [ElPalOJ. First note that the Stratonovich 
solution defined by (I28l)-(|29l)-(l3"0l) with the Wiener process (9l£/(f),3£/(0) 
coincides with the Ito solution because the vector fields s'm(q) d p and cos(g) d p 
commute. Finally, since cos 2 q + sin 2 q = 1, the process defined by dP = 
sin Q <MU (t) + cos Q d3 U (t) is distributed as the Wiener process in C( [0, 2k] , R) . 

□ 

Now we turn to the limit g/fm — > °°. In this limit, we know that the 
velocity components of the motions of jY particles also converge jointly in 
distribution to JY independent Wiener processes. The previous results then 
imply 

Theorem 4.3. Given JY initial data (q^Po) in R 2 (I < I < JY), such that 

1 — cos(<?q — qf) ) + c\pq — Pq I > pairwise for some c > 0, consider the 
resulting solutions to ^24\)-<[25\)-<[27\). Then given any K > 0, for g/fm — > 
°°, N — > °°, M — > oo, the jY -dimensional process p^ converges to an tri- 
dimensional Wiener process, and convergence is in law in C( [0, 2kK] , R^) 
with the topology of uniform convergence. 

This follows immediately from Theorem 3.1 in [ElPalO], using the brow- 
nian limit U and the continuous mapping theorem as in the proof of Propo- 
sition @3]just given. 

5. Interpretation of the results 

Our formulation of the limit theorem is rather formal, and our proof strat- 
egy differs from the more usual ones in the physics literature. 

This paper starts by reducing "noisy wave fields" to "white" ones in the 
dense spectrum limit, using a central limit theorem in the "wave field space" 
of functions U¥, as shown in sec.[3j Considering functions u and U is a way 
to get a handle on the limit process driving the particle motion, while it is 
harder to define directly the limit in terms of the noise "du/ dt". The dense 
spectrum limit is instrumental here to provide the many independent terms 
in the sum defining the wave field. 

The resulting wave field entails the brownian limit for the velocity of 
a single particle for < t < 2k [ElPalOJ. This single particle statement 
makes no reference to any velocity distribution function : we take a "tra- 
jectory" viewpoint on stochastic processes, and state a "diffusion process" 
limit rather than a "Fick equation" limit. 

The diffusion picture for JV particles also follows from our previous 
proof [ElPalO] that, if the wave field is a "periodic white noise", then par- 
ticles released in the resulting force field are independent in the srf — > °° 
limit. This independence between particles results from the fact that parti- 
cle velocities are a continuous martingale (viz., given the wave field history 
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and their own, their velocity increments have vanishing conditional expec- 
tation), from the fact that a martingale is completely characterized by its 
quadratic variation process (which eliminates the need for considering more 
than two particles jointly), and from the ergodicity of the random evolution 
of the relative velocity of any pair of particles. Estimates in HElPalOU are 
rather technical, and one may wish to revisit them to provide explicit rates 
of convergence. 

Our order of limits is important : first we take the dense spectrum limit 
s — > °o, then we let srf — > °°, and finally we consider jY > 1 and K > 1 , for a 
single realization of the wave field. Our convergence is in distribution with 
respect to the wave field random data. 

In contrast, usual arguments for the Fokker-Planck limit invoke a loss of 
memory for the particle motion, directly in terms of particle velocity. The 
gaussianity of the velocity distribution at a time t (given a Dirac at time 
0) is then seen as resulting from a central limit theorem with a sum over 
(rime-)successive independent increments. The quantity of interest (see e.g. 
eq. (9.32) in [HaWa04|) is often the (wave field) ensemble-averaged veloc- 
ity distribution function rather than the empirical distribution driven by a 
single wavefield. 

6. Perspectives 

We proved that the motion of jY particles in the field of random waves 
approaches a velocity-diffusion process in the dense spectrum limit. Our 
probabilistic proof highlights a central limit behaviour, while the Benisti- 
Escande proof stresses the elimination of correlations by appropriate changes 
of variables. In comparison with the latter proof, as well as with other 
derivations of the quasilinear limit, we show that uniformity of phases is un- 
necessarily strong an assumption : four-symmetry (viz. phase distribution 
invariant modulo tt/2) is sufficient. We also show that the wave amplitudes 
need satisfy only rather mild assumptions. 

Our proof uses the specific dispersion relation of Benisti and Escande, 
k m = k for all waves, and the regular spacing of phase velocities as o n = 
n/N. The first assumption enables the decomposition E m sin(g — (O m t + 
<p m ) = C(t) s'mq + S(t) cosq with coefficients C and S independent of q, and 
the second one permitted to use the large body of knowledge on random 
Fourier series. Relaxing these assumptions is physically desirable and will 
be considered in future work. 

In contrast with most earlier works in the plasma physics community, 
our formulation focuses on full particle trajectories, rather than one-particle 
distribution functions. In particular, the joint convergence theorem [473] sup- 
ports the familiar picture that the evolution of the empirical distribution 



16 



YVES ELSKENS 



,JY~ X Y&=\ «5( v — w/ m ) a - s - approaches the solution of the diffusion equa- 
tion d t f — d v Dd v f = for large JV : this law of large numbers, and fluc- 
tuations around it, require a further limit (Jf — > °°) to be discussed in the 
light of Ito's arguments [Ito83J. In substance, our Theorem 14.31 establishes 



for velocities what Lebowitz and Spohn [LeSp83 1 called Assumption C on 



the motion of particles in position space in order to derive Fick's law for 
self-diffusion. 

Another extension, in the case a„ = n /N, would be to allow K = kN with 
a fixed k in Theorem 14.3 [ for it would validate the diffusion picture for times 
beyond the discretization time Tdi sc '=2n/ (kAvy) = 2nN viz. the time scale 
after which the wave Fourier spectrum shows its discreteness. While some 
physical applications of the dense spectrum limit may be viewed as rejecting 
^disc t° infinity, the mathematical issue is interesting because of evidence 
that the diffusion description applies to the single-particle evolution over 
long times BeEs971lF^EIMlFJltl . 

Acknowledgements 

This work benefited from many discussions with D. Escande and mem- 
bers of equipe turbulence plasma, with E. Pardoux, and with participants 
to the 107th statistical mechanics conference at Rutgers. Stimulating com- 
ments by D. Benisti and D. Escande, and an explanation by A. Lejay are 
gratefully acknowledged, as are the careful reading and constructive com- 
ments by the anonymous referees. 

Appendix 

In Ref. [ElP alOII we introduced the auxiliary process (X t7 Y t ), describing 
the relative position and velocity of two particles evolving in the same wave 
field. This process solves 

(33) dX t = Y t dt , X = x, 

(34) dY t = sm(X t )dB t , Y =y, 

in the state space E = T x R\ {(0,0), (ic,0)}, where T = R/(2nZ) and B, 
is the standard brownian motion in C(R + ,R). We proved there in Proposi- 
tion 5.1 that, for any (x,y) E E, this process a.s. does not reach the points 
{(0,0), (7T,0)} in finite time. The proof in Ref. HElPalOI does not identify 
points modulo lit for their x component ; one can streamline it as follows. 

Proposition .1. For any (x,y) e E, M{t > : sin 2 (X t ) +Y t 2 = 0} = +°o 
a.s., andinf{0 > : limsup f ^ e - (sin 2 (X f ) +Y t 2 ) = +°°} = -f°° a.s. 

Proof Let R t = sin 2 (X,) + Y t 2 and define Z t = \ogR t . Denote by x either of 
these stopping times, corresponding respectively to Z t — > — °° and Z t — > +°°. 
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Then Ito calculus on [0, %\ yields 

(35) dsin 2 X f = (2sinZ r cosZ t )F r d/, 

(36) dF 2 = 2Y t sm(X t )dB t + sin 2 (X t ) at 

)S (X,) + sin 2 (* 



2Y t sin(X,) cos(X t ) + sin 2 (X ? ) Y 2 sin 2 (Z f ) 
dZ; = — — — -at — 2— ^ — -at 



(37 ) -f2^^dS,. 

Rt 

Noting that 2\ab\ <a 2 + b 2 and that | cos x\ < 1 yields the estimates Y 2 sin 2 (X,) < 
Rj/4 and 2F f sin(Z f )cos(X f ) + sin 2 (X r ) > -R t , so that on the time interval 
[0,T[ 



(38) 



3t r f 

Z t >Z --+ (p s dB s 
2 Jo 



where \<p s \ < 1. This ensures that Z t is bounded from below on any finite 
time interval since B t is bounded. Hence mf{t > : R t = 0} = +°° a.s. 
Similar upper estimates imply 

(39) Zt<Z + 2t+ f (p s dB s 



ensuring that Z t is bounded from above on any finite time interval. Hence 
inf{0 > : limsup f ^ e -i? f = +°°} = +°° a.s. □ 
The second claim of the present statement does not supersede Lemma 5.4 
of Ref. [ElPalO], which proves that Y t does a.s. not diverge as t — > °o. The 
present statement only proves that (X t , Y t ) remains in E for all t > a.s. 
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